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ABSTRACT 

The Galaxy appears to be richer in young, massive stellar clusters than previously 



known, due to advances in infrared surveys which have uncovered deeply embedded re- 
gions of star formation. Young, massive clusters can significantly impact the surround- 
ing interstellar medium (ISM) and hence radio observations can also be an important 
tracer of their activity. Several hundred cluster candidates are now known by examining 
survey data. Here we report on multiwavelength observations of six of these candidates 
in the Galaxy. We carried out 4.9 and 8.5 GHz VLA observations of the radio emission 
associated with these clusters to obtain the physical characteristics of the surrounding 
gas, including the Lyman continuum photon flux and ionized gas mass. Spitzer Infrared 
Array Camera observations were also made of these regions, and provide details on the 
stellar population as well as the dust continuum and polycyclic aromatic hydrocarbon 
emission. When compared to the known young, massive clusters in the Galaxy, the six 
cluster candidates have less powerful Lyman ionizing fluxes and ionize less of the H II 
mass in the surrounding ISM. Therefore, these cluster candidates appear to be more 
consistent with intermediate-mass clusters (10 3 -10 4 M ). 



INTRODUCTION 



Massive stars return a significant fraction of their mass to the interstellar medium (ISM) 
by means of stellar winds and supernovae. Massive stars are rare, have short lifetimes, and are 
predominantly observed in young, massive (< 20 Myr, > 10 4 M ) clusters. To date, few such 
clusters are known and well-studied in the Galaxy because many are located in the highly obscured 
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Galactic disk. A complete census of young, massive stellar clusters in the Galaxy is needed to 
determine their total number and distribution. Estima tes have shown that the Galaxy should host 
close to 100 of these clusters (jHanson &; Popescull2008l ). However, only a few such clusters are well 
known and studied. For example, the Arches, Quintuplet, and Central clusters in the Galactic 
center were uncovered with infrared spectro scopic observations , and are responsib le for heating 



and s haping much of the surrounding ISM (jCotera et al 



1996; 



1997, 



2001a, 



Figer et al.1 Il999al b: lLang et al 



2002j). In additio n, several other you ng, mas sive clusters in the Galaxy have been 



recen tly studied: RSGC 1 fc 2 (Davies et alJl2008f>. W49a (IHomeier fc Alyesl 120051 : Ide Pree et al 
19971). NGC 3603 dMelena et alJliood: borissova et alJbood : Ide Pree et al.lll999h . and Westerlund 



dBrandner et al.ll2008l : kothes k Doughertvll2007n . 



Over the last decade, surveys in the near-infrared and mid-infrared (2MASS and ffpiteer/GLIMPSE 



in pa rticular) have revealed as many as 1000 candidate clusters (e.g. iBica et al.ll2003l ; iMercer et al 



2005). A sample of 40 possible young, massive clusters were selected from these candidates using 



multiwavelength data, including radio and optical observations. Among these 40 sources, six were 
identified as havin g bright associated radio emission from the NRAO VLA Sky Survey (VLA-NVSS; 



Condon et al.lll998l ) . Here, we present a multiwavelength study of these six cluster candidates in the 



Galaxy. The radio observations in this study have higher sensitivity and resolution than the VLA- 
NVSS data, and provide greater insight into the properties of the gas and dust surrounding these 
clusters. The primary goal of the radio data is to determine the properties of the ISM surrounding 
the candidate massive clusters, and to compare them to the known young, massive clusters in the 
Galaxy. Estimates of the ionized mass and ionizing flux can be obtained with the radio data, from 
which the candidate sources can be categorized as supermassive or intermediate-mass clusters. The 
primary goal of the infrared data is to characterize the infrared emission surrounding the clusters 
as well as the infrared emission arising from cluster members. Comparison between the radio and 
infrared observations provides additional insight about the interstellar environment of the stellar 
clusters and the impact of the stars on the surrounding medium. 



2. CANDIDATE CLUSTERS: BACKGROUND INFORMATION 



Presented here is a summary of the background information for the six candidate massive 
stellar clusters. Table [T] lists the six candidate sources, including the most recent distance estimates 
and alternative identifiers. Note that the names used here have been shortened from names that 



originate in the Bic a & Dutra catalogs of infrared star clusters (jBica et al.ll2003l ; iDutra et al.ll2003l ; 



Dutra &: Bical l200ll ) . Sharpless (Sh2) n umbers are give n for the clusters that have been observed 



in the Sharpless catalog of H II regions (jSharpless!ll959l ). 
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2.1. BD52 



BD52 is a stellar cluster located in the H II region Sh2-187. iRusseil et al.l (|2007l ) used the 
spectrophotometric distance of a single star (B2.5V) to det ermine a distance of 1.44 ± 0.26 kpc to 
the stellar cluster. According to a multifrequency study bv lJoncas et al.l (|1992l ). the H II region is 



fairly young with a kinematical age of ~2xl0 5 yr. Optical and 1.4 GHz radio emission observed in 
this study suggest BD52 is enshrouded in its parental cloud, and that a BO zero-age-main-sequence 
star is the exciting source. The presence of a water maser demonstrated bv lZinchenko et al.l (|1998l ) 
provides evidence that this H II region could be undergoing star formation. 



2.2. BD65 



BD65 is a stellar cluster located in the H II region Sh2-209. iKlein et al.l (|2005l ) have studied 
this complex at millimeter wavelengths to look for evidence of embedded star formation. Their 
observations show that the region around BD65 can be described as a ring of cloud cores around a 
star cluster, and the complex may be a site of triggered star formation. IKlein et al.l (|2005l ) adopted 
a distance of 8.2 kp c, yet there is some discrepancy in t he distance to this source (estimates range 
from 4.9 to 9.8 kpc: lBica et aDbood : ICaplan et al.ll200oh . A distance to BD 65 of 9.8 kpc i s adop ted 
in this work based on the kinematic and photometric distance provided by C apian et al. (j200ol ). 



2.3. BD84 

This appears to be one of the first times BD84 has been thoroughly studied. However, it has 
been observed in several H II region surveys and has been associated with the H II re gion Sh2-283. 



Its dis tance has been determined spectrophotometrically to be 8.06 ± 0.30 kpc by IRusseil et al 



(120071 ). 



2.4. BD95 



The Galactic H II region Sh2-294, or BD95, has been studi ed in several pape rs, mostly H II 
region surveys, including a multiwavelength study of the region (jSamal et al.ll2007l ). The distance 
to BD95 has been estimate d several times; mos t recently it has been determined spectrophotometri- 
cally to be 3. 24 zfc 0.56 kpc (IRusseil et al.ll2007r). T he ionizing star is believed to be a star of spectral 
type ~B0V (jSamal et al.ll2007l ). lYun et al.l (J2008) confirm the ionizing source, and quote a spectral 
type of B0.5V. Their high resolution near infrared images reveal a double cluster morphology with 
one cluster fully embedded and optically invisible. 
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2.5. DB7 



DB7 is associated with an H II region known as Sh2-307. To the SW of DB7 is an ul- 
tracomp act H II region and photodissociation region (PDR) associated with the DB8 stellar 
cluster (IMessineo et al.l 120071 ). DB8's distance has been determined spectrophotometrically by 
Messineo et al.1 (|2007l ) to be 2.65 ± 0.4 kpc. We adopt this same distance for DB7. If the two 
clusters are asso ciated with the same H II complex, DB7 may have formed first and triggered the 
collapse of DB8 (|Messineo et alJhoOTl ). 



2.6. DBCL23 



The most recent detailed study of DBCL23 was done bv lComeron et al.l (|2005), who estimated 
an upper limit of 3.5 kpc for its distance. The spectrum of the nebula around this source indicates 
a temperature of 4 0,000 - 45,000 K (cor responding to an 05 - 07 star) for the ionizing stellar 
spectrum. However, IComeron et al.l (|2005l ) discovered a significant discrepancy between the ionizing 
flux supplied by a mid O-type star and that inferred from the optical and radio parameters. It is 
possible that DBCL23 is ionized by an external O-type star, although the ionizing source was not 
identified. Their results confirm that DBCL23 is not a massive star forming site. 



3. OBSERVATIONS AND IMAGE PROCESSING 

Radio observations were made at 8.5 and 4.9 GHz with the Very Large Array (VLA) telescope of 
the National Radio Astronomy Observatory (NRAO) H in August 2008 through observing program 
AL723. The VLA was in its D array configuration allowing for resolutions of 5-20". At both 
frequencies, J1331+305 was used as a flux density calibrator. Standard procedures for calibration, 
editing and imaging were carried out using the Astronomical Image Processing Software (AIPS) of 
the NRAO. Table 2 summarizes the image parameters for the different sources. All images were 
made with robust = weighting. 

Images of all six targets were taken in 3.6, 4.5, 5.8, and 8.0 /iin bands using the Infrared Array 
Camera (IRAC) onboard the Spitzer Space Telescope (Program ID 30734, PI Figer). The 5 position 
Gaussian dithering pattern was used, and the resulting images in each wavelength cover a 5' by 
5' area, with a pixel size of 1.2" . The total integration time was 471 s in each of the four bands. 
The Spitzer MOPEX tool was applied to the basic calibrated data to create mosaicked images of 
each band. An IDL adapted version of DAOPHOT was used to perform photometry on each band, 
and the results were cross-correlated. If a source was not present in one or more of the 4.5, 5.8, 
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or 8.0 11m. bands, we reported a null result for the specific band(s), rather than remove the source 
entirely. 

4. RESULTS 

In the following sections the morphology of the clusters and their surroundings in both radio 
and infrared images is discussed. The correlation between the observed radio and infrared emission 
in Figures Q3[7] is also examined. The wavelengths at which the observations were taken provide 
information about the type of emission produced by the stellar clusters. Ionized hydrogen in the 
clusters emits thermal free-free emission, which can be traced at radio frequencies. Emission from 
the thermal dust continuum is seen at infrared frequencies. The IRAC bands at 3.6, 5.8, and 8 
/im contain polycyclic aromatic hydrocarbon (PAH) features. PAH emission typically arises in 
PDRs and can be used to delineate the edge of the PDR. The 3.6 /im band is dominated by stellar 
emission, and the 8 /iin band contains primarily diffuse continuum and PAH emission. These two 
bands were used for the colorscale seen in Figures [T][7] with contours overlaid representing the radio 
emission at 8.5 and 4.9 GHz. Comparing the distribution of radio and infrared emission provides 
us with knowledge about the morphology and emission mechanisms of the stellar clusters and their 
H II regions. 

4.1. Radio and Infrared Images 

4.1.1. BD52 

Figure Q] shows BD52 with contours of radio emission (at 8.5 GHz (left) and 4.9 GHz (right)) 
overlaid on 3.6 and 8 /jm images. The radio emission associated with BD52 has an extended 
morphology. The 8.5 GHz radio contours in the left image show a curved, partial-shell structure 
of approximately 1' (thickness) x 3' (diameter) (0.4 pc x 1.3 pc at d = 1.44 kpc) in size, whereas 
the 4.9 GHz radio emission shows a complete shell-like structure as seen in the image on the right. 
Given the distribution of radio emission at 4.9 GHz, it appears that there must be missing flux at 

8.5 GHz. Indeed, the largest angular size of the VLA interferometer at 8.5 GHz in D-array is ~3'. 
This means that extended structures larger than this angular size (3') will not be detectable at this 
frequency. The 4.9 GHz emission associated with BD52 shows that the source has a diameter of 
about 4' (1.7 pc at d = 1.44 kpc), therefore, flux is missing in the 8.5 GHz image. Two components 
are visible in the 3.6 and 8 |im emission: the stellar sources and the diffuse emission that appears 
to be part of a larger infrared shell extending beyond the 4.9 GHz radio shell. The emission is 
more apparent at 8 ^m shown in the right image, and to the north looks more complete than at 

3.6 /im. The stellar cluster is visible in the 3.6 /im image on the left. An investigation of the stellar 
components of BD52 by Trombley et al. (in prep.) suggests a cluster mass of 1300 M and an age 
of 1 Myr. They conclude that several main sequence early B stars are producing the ionizing flux 
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which powers the H II region. 



In the SW corner of the image is a bright radio point source (at approximately RA, DEC 
(J2000): 01 23 00, 61 48 15). This radio point source has no counterpart in t he infrared and is 
known to be a bright non-thermal extragalactic radio source (ISnell Ballylll986l ). We calculate the 
spectral index of this compact source to be -2.6±0.07, which is in agreement with the non-thermal 
nature of the source. 

Figure [2] shows a small area from the center of BD52 (approximately 50" by 50", corresponding 
to 0.35 x 0.35 pc at d = 1.44 kpc) with 8.5 GHz radio contours overlaid on the 3.6 //m image. In 
the center of the stellar cluster, there is a point-like radio source that appears to be coincident 
with an infrared source. The 8.5 GHz flux is 0.3=L0.1 mJy. This association is suggestive and one 
interpretation is that the radio emission arises from a stellar wind in one or more of the members 
of the stellar cluster. Presuming the radio emission associated with the stellar source i s due to the 



massive stellar win d, it is possible to calc ulate the mass- loss rate using Equation (1) of lLang et al 
(|2001bl ). based on iPanagia &: Fellil ( 19751 ). Assuming a terminal wind velocity of ~1000 km s _1 
and a distance of 1.44 kpc, the mass-loss rate is approximately 0.35xl0 -5 M yr _1 . This mass-loss 



rate is somewhat higher than expected for an early B-type main sequence star ([Sternberg et al 



2003). This particular radio source also has counterparts in the near-IR (HST/NIC3), mid-IR 
(Spitzer/IRAC), and X-ray (CXO) (Trombley et al. (in prep)). 



4.1.2. BD65 

Figure [3] shows BD65 with contours of radio emission (at 8.5 GHz (left) and 4.9 GHz (right)) 
overlaid on 3.6 and 8 jim images. The 3.6 and 8 /im infrared emission associated with BD65 is 
primarily diffuse with denser cores of bright emission that form a shell-like structure around the 
central cluster located at RA, DEC (J2000): 04 11 08, 51 10 00. The shell of infrared emission 
is outlined by the 8.5 GHz radio emission with a diameter of ~2' (5.7 pc at d = 9.8 kpc) shown 
in contours in the left image. This emission also extends to the NE and to the SW of the central 
stellar cluster. The 4.9 GHz radio emission in the right image is more extended than the 8.5 GHz 
emission. The morphology defined by this radio emission differs from that seen in the 3.6 and 8 /im 
emission. The two brightest cores of infrared emission (e.g. RA, DEC (J2000): 04 11 05, 51 10 00; 
04 11 05, 51 10 30) are also noticeable in the radio emission, but are not resolved into two distinct 
components. There is another notable region south of the central cluster at RA, DEC (J2000): 04 
11 04, 51 08 00 that is visible at both radio and infrared frequencies. This compact region of bright 
emission approximately 0.5' (1.4 pc at d = 9.8 kpc) in size is not associated with the extended 
emission at 8.5 GHz, but is part of the extended emission at 4.9 GHz. 
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4-1.3. BD84 

Figure H] shows BD84 with contours of radio emission (at 8.5 GHz (left) and 4.9 GHz (right)) 
overlaid on 3.6 and 8 /um images. The radio emission at both frequencies for BD84 shows a 
symmetric distribution of contours around the central cluster about 1' (2.3 pc at d = 8.06 kpc) in 
size with surrounding diffuse emission. The emission at 3.6 and 8 //m shows a shell-like structure 
with bright edges, whereas the radio emission is more complete with a single central peak in 
emission. The 4.9 GHz radio emission in the image on the right contains a large amount of diffuse 
emission that extends well past the infrared emission with an approximate size of 3' x 2' (7.0 x 4.7 
pc at d = 8.06 kpc). The 8.5 GHz radio emission in the left image only extends NE of the cluster. 
The peak for the infrared emission lies west of the central cluster at about RA, DEC (J2000): 06 
38 26, 44 30. The peak in radio emission, however, appears closer to the center of the cluster at 
RA, DEC (J2000): 06 38 27, 44 45. The radio source at RA, DEC (J2000): 06 38 34, 45 00 
appears distinct in the 8.5 GHz radio emission, but is part of the extended emission at 4.9 GHz. 

4.1.4. BD95 

BD95 is shown in Figure [5] with contours of radio emission (at 8.5 GHz (left) and 4.9 GHz 
(right)) overlaid on 3.6 and 8 /xm images. The 4.9 GHz radio contours show a complete shell 
with a diameter of ~5' (4.7 pc at d = 3.24 kpc). There is corresponding 3.6 and 8 /xm infrared 
emission only along the eastern side where the brightest infrared emission can be seen. The 8.5 
GHz contours shown in the left image appear drastically different with emission only around the 
brightest infrared emission. As with BD52, it appears that flux is missing at 8.5 GHz, most likely 
due to the aforementioned largest angular size limitations of the VLA. The structure of the 3.6 and 
8 //m emission is complex in comparison to the simple shell-like 4.9 GHz emission. The infrared 
images consist of a bright region of emission at RA, DEC (J2000): 07 16 39, -09 25 45 with diffuse 
arm-like features emanating out to the NE and SW. This peak emission is host to an embedded 
cluster (see §2.4). The cluster containing the ionizing source is visible at 3.6 //m in the left image 
and is located NW of the cluster at RA, DEC (J2000): 07 16 33, -09 25 15. 

4.1.5. DB7 

The radio emission associated with DB7 is shown in contours (at 8.5 GHz (left) and 4.9 GHz 
(right)) overlaid on 3.6 and 8 /im images in Figure [6l The infrared emission in DB7 has a different 
morphology than the radio emission. The 8.5 and 4.9 GHz radio emission contains a denser circular 
region of contours around the cluster, and some surrounding diffuse emission. The densest radio 
contours which encompass the central cluster (seen at 3.6 /mi) have a diameter of ~2' (1.5 pc at d 
= 2.65 kpc). The 4.9 and 8.5 GHz radio emission have similar morphology, although the 4.9 GHz 
emission in the right image extends farther west. The 3.6 and 8 fim infrared morphology resembles 
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that of BD95 with areas of bright emission and diffuse arm-like structures. One of these structures 
which arcs out to the NW has corresponding radio emission, most noticeable at 8.5 GHz in the 
left image. Other than this structure, the radio and infrared are not well correlated. The emission 
peaks are offset with the infrared peak occurring at RA, DEC (J2000): 07 35 33, -18 45 30, and 
the radio peak occurring at RA, DEC (J2000): 07 35 35, -18 45 30. The associated source, DB8, 
can be seen in the right image at RA, DEC (J2000): 07 35 39, -18 48 45. 

4.1.6. DBCL23 

Figure [7] shows DBCL23 with contours of radio emission (at 8.5 GHz (left) and 4.9 GHz (right)) 
overlaid on 3.6 and 5.8 /im images. The radio contours at both frequencies are very similar, and 
both are well correlated to the infrared emission. There is only a slight difference in morphology 
between the infrared and radio with the radio emission being nlled-in, and the infrared emission 
more shell-like with a central cavity. The contours of 4.9 GHz emission shown in the image on the 
right give a diameter of approximately 2' (2 pc at d = 3.5 kpc) for DBCL23. The diffuse 5.8 /im 
emission seen in the right image extends farther than the 4.9 and 8.5 GHz radio contours. There 
are double emission peaks separated by about 30" visible in the radio and infrared emission. The 
location of these peaks in the radio emission is offset from the peaks in the infrared emission. 

4.2. Relationship Between Radio and Infrared Emission 

Examination of Figures [T][7] reveal some overall trends for the ISM surrounding the candidate 
clusters: (1) In all but two sources, the infrared emission (3.6 and 8 /im) extends beyond the 
associated radio emission. The two exceptions to this are BD84 and DB7, for which the 4.9 GHz 
emission extends farther than the infrared emission; (2) The morphology of the infrared emission for 
the cluster sources BD52, BD84, and DBCL23 is shell-like with a cavity and bright emission along 
the edges. The radio emission in BD52, BD84, and DBCL23 is also more symmetric and roughly 
follows the infrared emission. The shell-like morphology at both radio and infrared wavelengths 
allows for a more in-depth analysis to be carried out for these three sources (§5.1); (3) On the other 
hand, BD65, BD95, and DB7 are more complex and less symmetric in morphology at radio and 
infrared wavelengths; and (4) the peak intensity seen in the radio emission is offset from the peak 
in the infrared emission. This trend is also noticeable in the slice analysis performed on the three 
cluster candidates BD52, DBCL23, and BD84 (see §5.1). 
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4.3. Radio Continuum Properties of H II Regions 

4-3.1. Calculations 

We can estimate the physical properties of the radio continuum emission by examining the flux 
density of each source. Flux densities have been measured within the contour representing 3 times 
the rms level at 4.9 GHz as seen in Figure [HJ Approximate angular sizes have also been determined 
from the full width at zero intensity of slices taken across each source. The 3<r contour at 4.9 GHz 
was used to delineate the edge of the sources. The integrated flux density and angular size values 
are reported in Table El 

Based on the radio continuum parameters, physical properties of the ionized gas associated with 



the candidate clusters have been derived using the formulation of lMezger &: Henderson! (| 19671 ). In 
using these formulae, a uniform density, spherical, ionization bounded H II region with an electron 
temperature of 10,000 K and an abundance of ionized helium to hydrogen of 0.1 is assumed. Table 
H] presents the derived quantities from the radio continuum parameters shown in Table El Due to 
flux missing from a couple of the sources in the 8.5 GHz observations, an accurate measure of the 
flux density could not be obtained, thereby preventing a useful analysis at this frequency. Using the 
formulae, we derive a linear radius of the equivalent sphere for the H II region (r), electron density 
(n e ), emission measure (EM), total mass of ionized gas (M ), and the Lyman continuum photon 
flux (Nl„ c ). The Lyman continuum photon flux is the rate of photons emitted capable of ionizing 
the surrounding hydrogen. This value provides a measure of the cluster's power by quantifying its 
ability to impact the surrounding gas. 

The authors acknowledge that the sources presented in this study are not well described by a 



uniform density sphere. However, past studies (most notably the work done by lBalser et al.l ll995) 



have shown H II re gion properti e s bas ed on radio continuum emission to have relatively limited 



model dependence. iBalser et al.l (|1995l ) compared three models to describe the radio continuum 
structure of 11 Galactic H II regions. Physical parameters were derived from the models using differ- 
ent geometries (homogeneous sphere and a spherical Gaussian) as well as by using peak properties 
in the continuum intensity maps. Comparison between the three models showed no systematic 
differences in the derived electron density and emission measure, and resulted in an overall scatter 
of 15% between modeled parameters. In addition, the ionization properties which are of greatest 
interest in this analysis are a function of frequency, electron temperature, distance, and flux den- 
sity, all of which are model independent. The main challenge lies in accurately measuring the flux 



density. IBalser et al.l (|1995l ) discovered there could be differences in measured flux density of up to 
~50% for extended sources, which would result in about a 50% difference in derived ionizing flux. 
However, the motivation of this work is an order of magnitude comparison of ionizing photon fluxes 
between the sources in this study and known massive clusters in the Galaxy. Therefore, the results 
of this analysis are adequate for the primary goal of this study. It should be emphasized that the 
physical quantities reported in Table [H derived from the radio continuum parameters are meant to 
be used for relative comparison. 
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The last column in Table H] gives an estimate of the spectral type for each cluster based on our 



calcul ated Lyman continuum photon fluxes and model stellar atmospheres derived bv lSmith et al 
(|2002l ). No correction for extinction was made for these estimates. The stellar spectral types 
given are the types that would be assigned to a single ionizing source which produces the observed 
Lyman continuum flux. It is recognized that the ionizing photon rate calculated for the stellar 
clusters most likely arises from more than one cluster member, and not a single ionizing source of 
the estimated spectral type. However, the estimate allows for a comparison between the results 
of this analysis and those determined in previous studies. In general, there was close agreement 
between the estimated spectral types indicatin g these results are consistent with past observations 



(e.g. BD95; Samal et al. 



2007; 



Yun et al] 12008 ). 



4-3.2. Properties 

The physical parameters derived in Table H] show variation in physical characteristics of the H 
II regions associated with the cluster candidate sources. The equivalent radii for the H II regions 
range between ~1 and 9 pc, with the majority of sources having linear sizes of ~l-3 pc. The 
electron densities range from 10 cm -3 up to 140 cm -3 with an average of ~48 cm -3 . The EM 
differs for all sources ranging between ~10 3 pc cm" 6 and 6xl0 4 pc cm" 6 . The total mass of ionzed 
gas associated with the H II regions varies widely from 7 to 2800 M . The cluster source BD52 has 
the smallest ionized mass of 7.4 M Q . BD84, BD95, DB7, and DBCL23 have intermediate masses 
averaging around ~50 Mq, whereas BD65 contains a remarkable 2800 M of ionized gas. The 
Lyman continuum photon fluxes spread over two orders of magnitude from ~10 47 s -1 to 10 49 s _1 
with an average of ~6xl0 48 s _1 . This spread in ionizing photon rates corresponds to a spread in 
stellar spectral type from early B-type to mid O-type. The majority of the candidate sources have 
Lyman continuum photon fluxes that would be produced by single B0 type star. 

In this candidate cluster comparison, BD65 stands out as having greater values for most derived 
parameters. It clearly has the largest linear radius of 9.1 pc, much greater than the other clusters' 
average of ~2 pc. Similarly, BD65 has an H II mass that is nearly 70 times greater than the 
average of the other clusters' ionized mass. BD65 has the greatest Lyman continuum photon flux 
of 3.5xl0 49 s _1 , corresponding to an ionizing source of spectral type 03V. The physical paramters 
of the H II regions associated with these clusters depend on the distance. Therefore, the discrepancy 
in distance estimates to BD65 implies the physical properties derived would be less e xtreme if a 
smaller distance was adopted. For example, at a distance of 4.9 kpc ( Bica et al. 20031 ) BD65 has 



a linear radius of 4.5 pc, an ionized mass of 490 M Q , and a Lyman continuum photon flux of 
8.8xl0 48 s . Even at this smaller distance BD65 has the largest linear size, ionized mass, and 
ionizing photon flux of the candidate sources. Although small in size and mass as compared to 
BD65, DBCL23 has the largest electron density and emission measure. DBCL23 also has the second 
highest ionizing photon flux (2.7xl0 48 s" 1 ), which is about a magnitude less than that of BD65. 
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5. DISCUSSION 
5.1. Radio and Infrared Comparisons 



Watson et al.l (|2008l ) performed an analysis on infrared dust bubbles using radio and mid- 
infrared data for bubbles in the Galactic plane which led to a generalized picture of infrared bubbles 
as ionized gas with hot dust surrounded by a PDR containing interstellar gas, PAHs, and dust. 
The inner ionized gas is traced by radio free-free emission, and the thermal dust continuum can 
be traced with infrared emission, which often lies outside of the region of radio emission. The 
surrounding PDR is dominated by PAH emission features which can be seen in IRAC bands at 3.6, 
5.8, and 8 [im. Here, we carry out a similar analysis using our radio and infrared data for several 
of the sources that appear to be most symmetric: BD52, DBCL23, and BD84. 



5.1.1. BD52 

Figure [9] shows a slice through the BD52 shell at 4.9 GHz, 8.5 GHz (convolved to match the 
4.9 GHz) and 8 ^m. The location of this slice is shown by the blue line in Figure [TJ The slight 
difference in flux at 8.5 GHz compared to 4.9 GHz is most likely due to missing flux, as previously 
discussed in §4.1.1. Figure [9] shows that the radio emission peaks in two places near 200" and 
300" along the slice. The 8 fim emission exhibits the same double peaked structure, but peaks at 
~100" and 315" along the slice. This plot shows that the majority of the radio emission is contained 
within the 8 ^m emission indicating that the ionized gas lies primarily interior to the infrared shell 
(thermal dust). Figure [TUl shows the same slice across BD52, but for the 3.6, 4.5, and 5.8 /im data. 
The 3.6 and 4.5 fim bands are more sensitive to stellar emission, and show only faint continuum 
emission compared to the 5.8 /im emission. The 5.8 //m emission closely follows the structure of 
the 8 fim emission in Figure [9j The lower intensity emission between the two emission peaks at 5.8 
and 8 /im can be attributed to hot dust. The difference between this emission and the emission 
peaks at 5.8 and 8 /xm is likely due to the addition of PAH features in the peaks, which delineate 
the edge of the PDR around the cluster. 



5.1.2. DBCL23 

Figure HU shows a slice through the shell of DBCL23 at 4.9 GHz, 8.5 GHz (convolved to 4.9 
GHz) and 8 ^m. The location of this slice is shown in Figure A pair of peaks in the radio 
emission occur at ~90" and 120" along the slice. The 8 emission also shows double peaks that 
occur at ~65" and 130" along the slice. The majority of the radio emission is contained within the 
8 //m emission as it is in BD52 indicating that the ionized gas lies interior to the thermal dust. 
Figure [12] shows the same slice across DBCL23, but at 3.6, 4.5, and 5.8 ^m. The 3.6 and 4.5 fim 
emission roughly follows the 5.8 fim emission, although the source is significantly fainter at these 
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wavelengths. As in BD52, the 5.8 //m closely resembles the 8 fim emission. The dip in intensity 
between the emission peaks at these wavelengths is likely due to hot dust with the omission of 
PAHs. The PAHs are traced by the 5.8 and 8 //m emission in the two emission peaks. The spike 
at 120" seen in all three bands in Figure [T2l represents stellar emission. 



5.1.3. BD84 

Figure H3] shows a slice across the source BD84 at 4.9 GHz, 8.5 GHz (convolved to 4.9 GHz) 
and 8 /iin. The location of this slice is shown in Figure HI The radio emission for this source 
has a single peak near ~115-120" along the slice, and a FWHM of ~30". The 8 fim emission has 
peaks on either side of the radio emission, ~90" and 130". A possible interpretation is that the H 
II emission indicated by the peak is surrounded by a shell of PAHs and warm dust indicated by 
the two infrared emission peaks. As is evident in Figure [H the radio continuum emission is more 
extended than the infrared emission in BD84. This is only slightly apparent in the slices in Figure 
[13] where the 4.9 GHz emission is visible above the 8 fim emission just outside the two peaks. The 
different morphology of BD84 implies there is more ionized gas in this cluster system, and it is not 
all contained within the PDR. Figure [T4l shows the same slice across BD84 at 3.6, 4.5, and 5.8 /im. 
It can be seen in Figure [TH that the 5.8 /im band follows the 8 /xm emission very closely. The 3.6 
and 4.5 //m bands exhibit the same general structure, although the second peak at ~130" is much 
smaller, relatively. 



5.1.4- Slice Comparison 

The double peaked radio emission for the cluster candidates BD52 and DBCL23 implies a 
shell- like morphology for these H II regions. BD84 differs in this respect with a single peak of 
radio emission near the cluster surrounded by more diffuse extended emission. Lower intensity 
emission at 5.8 and 8 /im between two emission peaks for all three sources indicates PAHs are 
probably being destroyed by stellar radiation near the clusters. This also provides more evidence 
that the infrared emission peaks are due to PAH features, and represent the edges of the PDR. 
Our simplistic analy sis for these three sources is consistent with the infrared dust bubble model of 



Watson et al.1 (|2008l b 



5.2. Associated Stellar Sources: Infrared Color-Color Diagrams 



Photometry from Spitzer/IRAC can be used to i dentify objects wit h infrared excess, including 
evolutionary classes of young stellar objects (YSOs; lAllen et al.ll2004l ). YSOs serve as signposts 
for sites of ongoing star formation, and can be used to probe age differences in regions of extended 
star formation. In the simplest case, stars with infrared excess can be identified as candidate YSOs 
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(cYSOs), as the infrared excess is often indicative of a disk structure. Objects which are bright and 
extended in the 4.5 //m band, called extended green objects (EGOs), have a high as s ociati on rate 
with outflows from massive young stellar objects (MYSOs; Cvganowski et al. 20091 . boill ). Sites 



of massive star formation regions, e.g. in and around massive stellar clusters, can be selected by 
identifying a population of EGOs. 

Using simple assumptions, infrared photometry from Spitzer/IRAC can be used to construct 
color-color plots for sources in each field. Figures [T5ll20l show these plots for all six cluster candidates. 
In each plot, the left panel shows the -^3,6-4.5, M3.6-5.8 color-color plot. The right panel shows the 
distribution of different sources on the 8 fim image of each cluster candidate. Depending on the 
location of the sources on the color-color plot, they can be divided into three groups d esigned to 



probe different classes of objects: Group I, II and III. Following the procedure outlined in lZhu et al 



(|2009l ) we construct these groups as follows. Stars with little infrared color difference, e.g. main 
sequence stars, congregate around M3.6-4.5, -M3.6-5.8 = 0,0 (i.e., no excess color). Such objects are 
denoted as Group I in this work. A reddening vector corresponding to Ak = 10 is overplotted on 
each color-color plot. Fore- and background stars affected by differential reddening will follow the 
reddening vector, and are henceforth denoted as Group III. As stated earlier, excess emission at 
5.8 ^m can be indicative of strong PAH emission from circumstellar emission, e.g. a disk around a 
YSO, and are here referred to as Group II objects. In color-color space, these objects occupy the 
space defined by M3. 6-5.8 > 0.9 mag and -0.25 mag < M3, 6-4.5 < 0.25 mag. Stars with infrared 
colors in between Groups II and III (see, for example, the left panel of Figure PT6|) are also likely 
cYSOs. 

Active galactic nuclei can contaminate infrared color-color diagrams, though the contamination 
rate in this sample is likely low given the low latitude of these stellar clusters. Here we discuss the 
infrared excesses seen toward each of the target sources. Only objects with valid photometry in 
all four bands are considered. Examination of the distribution of Group II sources indicates that 
all six clusters are sites of ongoing star formation. Younger, more deeply embedded protostars are 
undetectable at these wavelengths. 



5.2.1. BD52 

The IRAC color-color diagram of BD52 is shown in the left panel of Figure [15j The spatial 
distribution of Group II and Group III stars is displayed in the right panel. Group III stars are 
fairly evenly distributed across the face of the cluster, though a dearth of these objects are seen in 
the direction of bright 8.0 /jm emission peaks. As these sources are distributed along the reddening 
vector, it is likely that their colors are affected by varying extinction. The 8.0 fj,m emission indicates 
regions of high local extinction, making it difficult to observe background stars in these areas. Only 
one object appears in the Group II range of colors. This object (shown as a diamond in the right 
panel of Figure fl~5"j) is located at the edge of the shell-like 8 /mi emission. 
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5.2.2. BD65 

Figure [16] displays the IRAC color-color diagram of BD65, left, along with the spatial dis- 
tribution of objects with infrared excess, right. Group II objects, indicated by a diamond in the 
right panel of Figure [16] appear to be distributed preferentially along ridges of the bright, extended 
8.0 //m emission surrounding BD65. The few Group III objects in this region are associated with 
point sources in the 8.0 ^m image, suggesting that the majority are foreground stars distributed 
along the line of sight. Stars farther from the reddening vector, falling mostly into Group II, are 
indistinguishable from cYSOs. 

5.2.3. BD84 

The left panel of Figure [T7] shows the color-color diagram for the sources detected, while the 
right panel shows the distribution of objects overlaid on the 8 /iin image. The number of infrared 
sources near BD84 is much lower than for the other target regions. There are only ~35 sources 
with valid photometry in all four bands; ~100 with valid photometry in 3 bands. The majority of 
sources in the field of view are likely to be fore- and background stars. Group III objects follow 
the reddening vector closely; about half of these sources are located preferentially along a ridge of 
infrared emission to the West of BD84 and at further distances from the bright infrared emission. 
Group II sources are more tightly clustered around BD84 and the surrounding bright infrared 
emission. It is likely that many sources associated with the cluster are undetected. 

5.2.4. BD95 

The left panel of Figure [18] shows the color-color diagram for the sources near BD95 and the 
right panel shows the distribution of the Group II and III sources overlaid on the 8 fj,m emission. 
Both Group II and Group III stars are preferentially located close to the stellar cluster BD95 and 
along ridges of infrared emission. The majority of these objects are likely cYSOs, highlighting sites 
of ongoing star formation in and around BD95. 

5.2.5. DB7 

Group II sources in DB7 are located along filaments and knots of bright 8.0 /xm emission, 
which can be seen in the right panel of Figure [19] These stars are also located to the SE, in the 
direction of the stellar cluster DB8. The majority of Group III sources are located fairly far from 
the reddening vector, and seem to be concentrated in and around the two clusters, making it likely 
these are cYSOs. 
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5.2.6. DBCL23 

Inspection of the right panel of Figure 1201 shows that Group II sources are located preferentially 
towards the bright, diffuse 8.0 fim emission surrounding DBCL23. The Group III objects located 
outside of DBCL23 are likely reddened background sources, while those within the cluster are 
associated with bright point sources, making these sources likely cluster members. 



5.2.7. Overview 



Inspection of the spatial distribution of objects with infrared excess in each of the six clusters 
presented in this study yield similar results. Each cluster contains candidate YSOs in regions of 
extended 8.0 /im emission. It is possible that the distribution of infrared excess sources results 
from star formation triggered by the expansion of bubble-like H II regions, but it can also be 
argued that already formed protostellar objects are being exposed by the evolution of the central 
systems. A more robust analysis of the YSO population of each cluster would be possible at 
wavelengths between those investigated in this study, e.g. Herschel and/or ALM A observations 



Thes e observations would also prove useful as probes of triggered star formation (see lZavagno et al 



2010 for one such study) . Maser activity is often assoc iated with shocked emission such as outflows 



from young stellar objects (e.g. lAnglada et al.lll996l ): in the case of BD52 the possible EGO is 



cut off by the edge of the IRAC f rame, so a firm association can not be made. The maser- YSO 
association is also pointed out by ICyganowski et al.1 1)20081 . l2009i ) who focus on the association 
between EGOs and masers. The Spitzer/IRAC observations of each cluster were inspected by eye 
for EGOs, yielding no obvious sources. A lack of detection of MYSOs (bright, extended sources 
in the 4.5 /im discussed above) is consistent with the general picture for all six clusters: young, 
partially embedded, intermediate mass stellar clusters. The lack of MYSO detection also implies 
the candidate clu sters are not sit e s of o ngoing massive star formation. For DBCL23, this result 
is consistent with IComeron et al.1 (|2005l ) 's confirmation that it is not a massive star forming site 
(§2.6). 



5.3. Comparison to Known Young, Massive Clusters in the Galaxy 

Table [5] lists physical properties of 14 young, massive clusters in the Galaxy that have been 
well-studied with which the properties of the candidate clusters can be compared. The properties 
most relevant to this study are the H II mass and the Lyman continuum photon flux, because they 
quantify the impact of these stellar clusters on the surrounding ISM. It is apparent that in general 
the average mass of ionized gas implied by the radio observations for the candidate sources is below 
the well-known massive clusters. With the exception of BD65, all other cluster candidate sources 
have ionized masses of 7-70 M©. BD84 and DBCL23 are at the upper end of this range with ionized 
masses of 70 and 56 M , which are comparable to the ionized mass associated with Westerlund 
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1 (jKothes &: Doughertvl 120071 ). BD65 stands out as having a very large ionize d mass of 2800 M , 
This mass is comparab le to that of the W49a and Quintuplet H II regions (jde Pree et ajj [1993 ; 
Lang et alJll997L boOlj ). 



The cluster candidate sources have ~2-3 orders of magnitude smaller Lyman continuum photon 
fluxes when compared to these other well-studied massive Galactic clusters, which have a typical 
value > 10 50 s . Our largest value is 3.5xl0 49 s _1 for BD65, and even this val ue is almost an 



order of magnitude less than the next closest value of ~10 50 s 1 for NGC 3603 (jde Pree et al 



1999). The candidate cluster DBCL23 is also on the higher end of our Ni yc range with a value of 
2.7x l0 48 s -1 . This i s comparable to DBS2003 45 which had a Lyman continuum flux of 8.2 xlO 48 
s" 1 (|zhu et aliboosh . 



The impact on the ISM from massive clusters such a s Arches and Quintuplet results in sizable 
structures 2 pc and 10 pc in size (|Lang et al.lll997ll2001al ). The size of the H II regions around the 
cluster sources BD65 and BD84 fall into this range. However, the linear radii of the H II regions 
around the rest of the candidates are more analogo us to the sizes of t he actual stellar clusters for 
some of the comparison clusters (e.g. RSGC 1 & 2; lDavies et al.l 1200a ). 

Past studies have predicted the Galaxy to harbor about 1 00 massive (~10 M ) stellar clusters 
similar to the comparison clusters (IHanson fc Popescull2008l ). However, using the evidence stated 
above, we conclude that the candidate clusters do not re present the most mass ive stellar clusters 
in the Galaxy, but are rather intermediate mass clusters. iMessineo et al .1 (|2009i ) extrapolated the 
number of known clusters with a mass larger than 10 4 M to conclude that there are about 200 
clusters with masses > 10 3 M . The stellar clusters presented in this study may belong to this 
group of intermediate mass clusters in the Galaxy. 



6. CONCLUSIONS 



Observations have been carried out using Spitzer/IRAC at 3.6, 4.5, 5.8 and 8.0 /jm and VLA 
radio observations at 8.5 and 4.9 GHz to determine the properties of six candidate stellar clusters 
and their surrounding H II regions. The following conclusions have been drawn from this study: 



1. Image overlays of all the sources were created to compare the distribution of infrared and 
radio emission. The 3.6 and 8 (im emission was found to be more diffuse than the radio 
emission in the majority of the cluster sources. In addition, the peak intensity in the radio 
emission is often offset from the peak intensity in the infrared emission. 

2. Radio continuum properties of the H II regions led to an estimate of the sources' ionized 
mass and Lyman continuum photon flux. BD65 is the most massive of the candidate clusters 
having an ionized mass of ~2800 M and ~10 49 . DBCL23 was also found to have a 
fairly large ionizing flux of ~10 48 . The physical parameters of the sources from the radio 
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continuum were compared to H II parameters of other well-studied young, massive Galactic 
clusters. 

3. A more complete multiwavelength analysis of the morphology in BD52, DBCL23, and BD84 
was carried out. The emission mechanisms in the clusters and surrounding ISM could be 
explored by investigating the intensity of emission at different wavelengths across a slice 
through the sources. A general picture emerged from this analysis: ionized gas surrounded 
by bright warm dust and PA H emission. This is consistent with other studies of infrared dust 



bubbles (jWatson et aU 12008 ) 



4. Infrared photometry of stellar members of the candidate clusters was used to search for can- 
didate YSOs. The existence of YSOs provides evidence of ongoing star formation. Candidate 
YSOs were found in all the cluster candidates using color-color plots. A lack of detection of 
MYSOs is consistent with the general picture for all six clusters: young, partially embedded, 
intermediate mass stellar clusters. 

5. Comparison with other young, massive clusters (e.g. Arches, Quintuplet, Westerlund 1) leads 
to the conclusion that the stellar cluster candidates presented in this study are not massive 
( > 10 4 Mq) clusters. More likely, they represent examples of Galactic intermediate mass ( 
> 10 3 M ) stellar clusters. 
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Fig. 1. — BD52: The left image shows radio contours at 8.5 GHz corresponding to levels of 3, 6, 
9, 15, 30, 60, 120, 200, 300, 450 times rms noise of 3.65xl0 -5 Jy beam -1 overlaid on color scale 
representing Spitzer IRAC data at 3.6 /im. The right image shows radio countours at 4.9 GHz 
corresponding to levels of 3, 5, 7, 9, 12, 15, 18, 50, 100, 200 times rms noise of 3.13xl0 -4 Jy 
beam -1 overlaid on color scale representing Spitzer IRAC data at 8 /im. The blue lines represent 
the slice taken through the shell seen in Figures l9l and [TUl The synthesized beam of the radio image 
is shown in the lower left corner of each plot. 
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Fig. 2. — BD52 center with radio contours at 8.5 GHz corresponding to levels of 3, 6, 9, 12, 15 
times rms noise of 3.65 xlO -5 Jy beam -1 overlaid on color scale representing Spitzer IRAC data at 
3.6 /mi. The synthesized beam of the radio image is shown in the lower left corner of the plot. 
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Fig. 3. — BD65: The left image shows radio contours at 8.5 GHz corresponding to levels of 3, 6, 9, 
15, 25, 35, 45, 60, 75, 90, 120, 150, 180, 220, 260 times rms noise of 1.56xl0~ 4 Jy beam" 1 overlaid 
on color scale representing Spitzer IRAC data at 3.6 /im. The right image shows radio contours at 
4.9 GHz corresponding to levels of 3, 8, 15, 22, 32, 45, 60, 80, 100, 150 times rms noise of 5.52 xl0~ 4 
Jy beam -1 overlaid on color scale representing Spitzer IRAC data at 8 ^m. The synthesized beam 
of the radio image is shown in the lower left corner of each plot. 
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Fig. 4. — BD84: The left image shows radio contours at 8.5 GHz corresponding to levels of 3, 6, 
9, 15, 25, 35, 50, 70 times rms noise of 4.25 xlO -5 Jy beam -1 overlaid on color scale representing 
Spitzer IRAC data at 3.6 fim. The right image shows radio contours at 4.9 GHz corresponding to 
levels of 3, 6, 9, 15, 25, 35, 50, 70 times rms noise of 6.75xl0 -5 Jy beam -1 overlaid on color scale 
representing Spitzer IRAC data at 8 //m. The blue lines represent the slice taken through the shell 
seen in Figures [13] and [141 The synthesized beam of the radio image is shown in the lower left 
corner of each plot. 
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Fig. 5. — BD95: The left image shows radio contours at 8.5 GHz corresponding to levels of 3, 5, 7, 
9, 14, 18, 22 times rms noise of 5.91 xlO -5 Jy beam -1 overlaid on color scale representing Spitzer 
IRAC data at 3.6 //m. The right image shows radio contours at 4.9 GHz corresponding to levels of 
3, 5, 7, 9, 12, 16, 20 times rms noise of 2.81 xlO -4 Jy beam -1 overlaid on color scale representing 
Spitzer IRAC data at 8 /im. The synthesized beam of the radio image is shown in the lower left 
corner of each plot. 
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Fig. 6. — DB7: The left image shows radio contours at 8.5 GHz corresponding to levels of 3, 6, 
9, 15, 30, 45, 60, 90, 120, 160 times rms noise of 1.06xl0 -4 Jy beam -1 overlaid on color scale 
representing Spitzer IRAC data at 3.6 fim. The right image shows radio contours at 4.9 GHz 
corresponding to levels of 3, 6, 9, 15, 25, 40, 60, 100, 150, 200, 250 times rms noise of 1.94xl0 -4 
Jy beam -1 overlaid on color scale representing Spitzer IRAC data at 8 /im. The synthesized beam 
of the radio image is shown in the lower left corner of each plot. 
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Fig. 7. — DBCL23: The left image shows radio contours at 8.5 GHz corresponding to levels of 3, 8, 
18, 32, 50, 80, 120, 145, 170, 225, 290, 325, 375 times rms noise of 1.34xl0" 4 Jy beam" 1 overlaid 
on color scale representing Spitzer IRAC data at 3.6 /im. The right image shows radio contours 
at 4.9 GHz corresponding to levels of 3, 6, 15, 30, 45, 60, 90, 120, 150, 200, 250, 300, 350, 400 
times rms noise of 2.95xl0 -4 Jy beam -1 overlaid on color scale representing Spitzer IRAC data at 
5.8 fim. The blue lines represent the slice taken through the shell seen in Figures [IT] and [12j The 
synthesized beam of the radio image is shown in the lower left corner of each plot. 
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Fig. 8. — Images show the 4.9 GHz radio continuum with contours at 3, 9, 18, 32, 80 times rms 
noise in each source. Flux densities were integrated for each source using the 3a contour. A 4' bar 
with the corresponding linear size at the source's distance has been included to demonstrate the 
relative size scale. 
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Fig. 9.— BD52 slice near DEC = 61°53' at 8 fxm (solid line), 4.9 GHz (dotted line), and 8.5 GHz 
(dashed line). The position of the slice is shown exactly in Figured! The 8 fim intensity has been 
scaled down by a factor of 10 5 . 
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Fig. 10.— BD52 slice near DEC = 61°53' at 3.6 (solid line), 4.5 (dotted line), and 5.8 fjm (dashed 
line). The position of the slice is shown exactly in Figured! 
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Fig. 11— DBCL23 slice near DEC = 44°15'15" at 8 fxm (solid line), 4.9 GHz (dotted line), and 
8.5 GHz (dashed line). The position of the slice is shown exactly in Figure [71 The 8 /jm intensity 
has been scaled down by a factor of 5xl0 3 . 
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Fig. 12.— DBCL23 slice near DEC = 44°15 / 15" at 3.6 (solid line), 4.5 (dotted line), and 5.8 fj,m 
(dashed line). The position of the slice is shown exactly in Figure [71 
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Fig. 13.— BD84 slice at DEC = 0°44'33" at 8 /an (solid line), 4.9 GHz (dotted line), and 8.5 GHz 
(dashed line). The position of the slice is shown exactly in Figured! The 8 fim intensity has been 
scaled down by a factor of 2xl0 4 . 
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Fig. 14.— BD84 slice at DEC = 0°44'33" at 3.6 (solid line), 4.5 (dotted line), and 5.8 (dashed 
line). The position of the slice is shown exactly in Figured! 
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Fig. 15. — Left: IRAC color-color diagram for BD52. Right: locations of stars showing 5.8 
excess (Group II) and strong foreground extinction (Group III), indicated with diamonds and 
squares, respectively. 




Fig. 16. — Left: IRAC color-color diagram for BD65. Right: locations of stars showing 5.8 ^m 
excess (Group II) and strong foreground extinction (Group III), indicated with diamonds and 
squares, respectively. 
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Fig. 17. — Left: IRAC color-color diagram for BD84. Right: locations of stars showing 5.8 /im. 
excess (Group II) and strong foreground extinction (Group III), indicated with diamonds and 
squares, respectively. 




Fig. 18. — Left: IRAC color-color diagram for BD95. Right: locations of stars showing 5.8 ^m 
excess (Group II) and strong foreground extinction (Group III), indicated with diamonds and 
squares, respectively. 
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Fig. 19. — Left: IRAC color-color diagram for DB7. Right: locations of stars showing 5.8 /iin excess 
(Group II) and strong foreground extinction (Group III), indicated with diamonds and squares, 
respectively. 




Fig. 20. — Left: IRAC color-color diagram for DBCL23. Right: locations of stars showing 5.8 
/im excess (Group II) and strong foreground extinction (Group III), indicated with diamonds and 
squares, respectively. 
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Table 1. Candidate Clusters 



Source 


RA,DEC 


Distance 


Alternative 


Name 


(J2000) 


(kpc) 


Identifiers 


BD52 


01 23 06.0, 61 51 32.0 


1.44±0.26 a 


Sh2-187, LBN 126.70-00.80, IRAS 01195+6136 


BD65 


04 11 09.4, 51 10 01.0 


9.8 b 


Sh2-209, G151. 59-0.23, IRAS 04073+5102 


BD84 


06 38 27.8, 00 44 39.0 


8.06±0.30 a 


Sh2-283 


BD95 


07 16 33.0, -09 25 26.0 


3.24±0.56 a 


Sh2-294, LBN 1032 


DB7 


07 35 33.7, -18 45 39.0 


2.65±0.4 C 


Sh2-307, ESO560-1, SF049 (DB8) 


DBCL23 


20 45 37.4, 44 15 14.0 


<3.5 d 


G84.0+0.8, ECX6-38, 3C 423 



HRusseil et al. I d2007tl 
Kaplan et al] fcoool l 
* lMessineo et al ] l|2007h 
^Comeron et alj l|2005n 



Table 2. Table for Radio Images 



Source 


Frequency 


Resolution 


rms 


Name 


(GHz) 


(arcsec) 


(mjy beam -1 ) 


BD52 


4.9 


17.93 x 11.61 


0.313 


BD52 


8.5 


8.77 x 6.86 


0.0365 


BD65 


4.9 


16.60 x 11.41 


0.552 


BD65 


8.5 


10.75 x 6.88 


0.156 


BD84 


1.9 


29.15 x 11.78 


0.0675 


BD84 


8.5 


14.68 x 7.37 


0.0425 


BD95 


4.9 


19.24 x 12.76 


0.281 


BD95 


8.5 


12.14 x 6.97 


0.0591 


DB7 


4.9 


24.27 x 11.77 


0.194 


DB7 


8.5 


15.66 x 6.93 


0.106 


DBCL23 


4.9 


24.78 x 12.30 


0.295 


DBCL23 


8.5 


15.70 x 6.98 


0.134 
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Table 3. Radio Continuum Parameters 



Source 


Flux Density at 4.9 GHz 


Angular Size a 


Distance 


Name 


(mjy) 


(arcsec) 


(kpc) 


BD52 


503±25 


250 x 250 


1.44±0.26 b 


BD65 


4100±200 


350 x 200 


9.8 C 


BD84 


42±3 


175 x 120 


8.06±0.30 b 


BD95 


279±20 


220 x 220 


3.24±0.56 b 


DB7 


511±30 


200 x 175 


2.65±0.4 d 


DBCL23 


2440±100 


120 x 140 


<3.5° 



a Approximate angular sizes determined using the full width at zero 
intensity of slices taken across each source delineated by the 4.9 GHz 
3<r contours. 

^Russeil et ah | (120071') 

fcaplan et all feOOCh 

^Messineo et al.l J2007h 

jComerdn et al.l (|2005l') 



Table 4. Physical Quantities of H II Regions Derived from Radio Continuum 



Source 


Linear Radius 


Electron 


Emisson 


Ionized 


Lyman 


Equivalent 


Name 


of Sphere 


Density 


Measure 


Mass 


Photon Flux 


Spectral 




(pc) 


(cm -3 ) 


(pc cm -6 ) 


(M ) 


(photons s _1 ) 


Classification 81 


BD52 


1.3 


37 


2.7X10 3 


7 


9.3X10 46 


B0.5V 


BD65 


9.1 


37 


2.6X10 4 


2800 


3.5X10 49 


03V 


BD84 


1.1 


10 


8.7xl0 2 


70 


2.4xl0 47 


B0V 


BD95 


2.5 


22 


2.5xl0 3 


35 


2.6xl0 47 


B0V 


DB7 


1.7 


1.3 


6.4xl0 3 


22 


3.2xl0 47 


B0V 


DBCL23 


1.6 


110 


6.4xl0 4 


56 


2.7xl0 48 


08V 



Note. — Quantities calculated assuming T e = 10,000 K and Y + = 0.1. 

a Spectral type estima t es of exciting stars are derived from Nr,ye using model stellar atmosphere 
results from lSmith et all hooj ) 



Table 5. Cluster Comparison Parameters 



Source 
Name 


Distance 


Galactic Coordinates 
I b 


(Mvr) 


Cluster Radius 
(dc) 


Cluster JVIeiss 

noooxMrO 


H II Mass 


Nr 

1 ' Lyc 

is' 1 ) 


Rcfcrei i< os 


Westerlund 1 


3.55 


339.55 


-0.4 


3-5 


1.0 


45 


47-53 




1.2 


Westerlund 2 


2.8 


284.27 


-0.34 


1-3 




7 




6xl0 50 


3,4 


RCW49 


8 


284.31 


-0.33 








... 


11.5x10 s0 


3 


RSGC1 


6.6 


25.27 


-0.16 


10-14 


1.5 


30 






5 


RSGC2 


5.8 


26.19 


-0.07 


14-20 


3.2 


40 






5 


W49a 


11.4 


43.17 


0.0018 


0.3-2 




50-70 


> 238 


2.4xl0 50 


6,7 


Central 


8 


0.00 


0.01 


3-7 


0.23 


20 




3.2 x 10 s0 


8,9 


Arches 


8 


0.12 


0.02 


2-3 


0.23 


< 70 


5480 


4x10 s1 


10,11 


Quintuplet 


8 


0.16 


-0.06 


3-5 


1.0 


20 


240 


7.9x10 s0 


9,12 


NGC 3603 


7.6 


291.63 


-0.53 


2.5 


0.7-1.5 


10-16 


10000 


1x10 s0 


13,14 


CL 1806-20 


8.7 


10.00 


-0.24 


3-4.5 




3 






15,16 


Glimpse 30 


7.2 


298.76 


-0.41 


4-5 


1.36 


3 






8 


DBS2003 179 


7.9 


347.58 


0.19 


2-5 


0.2 


7 






17 


DBS2003 45 


4.5 


283.88 


-0.91 


5-8 




1 




8.2xl0 48 


18 
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